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Based on a recent theory (Coppi, Nuclear Fusion, 42, 1, 2002) of spontaneous 
toroidal rotation in tokamaks (Lee et al, Phys Rev Lett, 91, 205003, 2003) and in 
astrophysical accretion disks, we propose that an analogous process could be at play 
also in the Earth space environment. We use fully kinetic PIC simulations to study 
the evolution of drift instabilities and we show that indeed a macroscopic velocity 
shear is generated spontaneously in the plasma. As in tokamaks, the microscopic 
fluctuations remain limited to the edge of the plasma channel but the momentum 
spreads over the whole macroscopic system. 



I. INTRODUCTION 

The understanding of the interplay between microinstabilities and macroscopic system- 
scale processes poses a formidable challenge to current plasma physics research. The present 
paper attempts to attract the attention of the plasma physics research community to a possi- 
ble analogy between space plasmas and tokamaks in the hope that a fruitful interdisciplinary 
exchange can lead to a new a deeper understanding. 

Among the small scale kinetic processes, the drift instabilities driven by density gradients 
have attracted considerable attention [l|. It has long been known that the small scale electric 
field and magnetic field fluctuations induced by such waves can act as an effective collision 
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by altering, the momentum of plasma particles in a manner analogous to regular classical 
collisions 2]. Such processes have been termed anomalous in the sense that when modeled 
within the framework of classical transport or fluid frameworks appear to lead to a much 
higher collisionality than is actually present 3|]. But such an approach focuses on the small 
scale momentum exchanges and sees the microinstabilities as acting in a localized fashion 
akin to collisions. Instead, a number of recent findings suggest that drift instabilities can 
act to alter momentum on a global macroscopic scale. 

Recent tokamak research has incurred in a new crucial experimental finding: the sponta- 
neous creation of toroidal rotation. In a number of experiments (including Alcator C-Mod, 
DIII-D, TCV, Tore-Supra, JET, TFTR, COMPASS) the plasma has been observed to ac- 
quire a toroidal rotation profile spontaneously in shots where no external input of toroidal 
momentum was imposed (see ref. N, \& and references therein). Researchers at Alcator 
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C-Mod have investigated the issue in depth [4], 17|, |8J. We refer the reader to their work for 
details and we limit the present discussion to highlight the findings most relevant to the 
present work. First, the generation of toroidal rotation is spontaneous and no external in- 
jection of angular momentum could be identified. Second, the source of angular momentum 
appears to be located in the edge zone where angular momentum is observed to form first. 
Later during the shot, toroidal angular momentum is transported inward at a pace faster 
than predicted by neoclassical diffusivity. 

To explain the observed spontaneous angular rotation in tokamaks, a number of possi- 
bilities have been proposed. We refer the reader again to the existing literature from the 
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Alcator C-Mod team for a discussion of all possibilities [8( and we limit the scope here to the 
particularly promising theoretical model proposed by Coppi [9]. Coppi noted an intriguing 
analogy between the spontaneous rotation observed in tokamak experiments and angular 
momentum exchange in accretion disks in astrophysical systems (such as forming stars or 
accreting black holes). The same core mechanism is at play in both systems: the exchange 
of angular momentum from the action of microscopic drift instabilities caused by the pres- 
ence of pressure gradients. When applied to tokamak plasmas, Coppi's model identifies the 
source of angular momentum at the edge of the plasma column where fluctuations due to 
edge-localized pressure- driven modes scatter particles out of confinement and transfer angu- 
lar momentum to the wall, thereby producing a net momentum in the plasma by reaction. 
The momentum created at the plasma edge is then transported inward and leads to the 
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observed profiles of toroidal angular velocity measured in Alcator C-Mod. 

In the present work we revisit the model proposed by Coppi and extend the analogy to 
a different class of plasmas: the space plasmas present in the near Earth environment. The 
usefulness of this analogy is that space plasmas are amenable to direct satellite observation 
of small scale processes and a wealth of data is available from past, current and planned 
missions. The paper is organized as follows. In section 2, we summarize the processes 
leading to the spontaneous formation of velocity shear in space plasmas. In Section 3 we 
present our simulation approach to the problem and in Section 4 we present the results of 
our simulation study and compare them with satellite observations. Conclusions are drawn 
in Section 5. 



II. LARGE SCALE VELOCITY SHEAR FORMATION IN SPACE PLASMAS 



The spontaneous formation of large scale velocity profiles has been recently discovered 
in space plasmas near the Earth jioillli 12] . The Earth is surrounded by a magnetic field 



that in its interaction with the solar wind forms a complex and dynamical structure, called 
magnetosphere (l3j. In the direction opposite to the Sun, approximately on the equatorial 
plane, the magnetic field is stretched to form the magnetotail where a current sheet is present 
and the magnetic field reverses sign. The northern part of the current sheet is magnetically 
linked to the north pole while the southern part is linked to the South pole. The current 
sheet covers a wide angular span in the night side of the Earth. Using the Geocentric Solar 
Magnetospheric (GSM) coordinates, the toroidal direction moving around the Earth in the 
magnetic equatorial plane is called the dawn-dusk direction, the vertical direction along the 
magnetic poles is the north-south direction and the third direction is from the Sun towards 
the Earth. 

Recent studies of the current layer in the magnetotail have uncovered a previously unex- 



pected result 
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ll| : microinstabilities on the scale of the electron gyroscale (of the order of 
the Km) developing in the edge of the current layer exchange momentum with the plasma 
in the dawn-dusk direction not only in microscopic regions, but extending over the whole 
current layer, spanning distances that equal tens of thousands of Km. 

The existence of density gradient driven instabilities in the lower-hybrid range (lower- 
hybrid drift instability, LHDI) has been an established fact for a long time, confirmed by 



numerous missions dating back from the early times of space exploration [14|. But it was 



always believed that such instabilities would affect the overa 



1 system in the form of sub- 



stitute collisions leading to anomalous transport effects The fact is, instead, that 

they lead directly to plasma-waves exchanges of momentum on the macroscopic scale of the 
whole current layer. 

The analogy with the tokamak experiments summarized above is obvious: in both cases 
the toroidal (or dawn-dusk) flow is generated spontaneously in the edge of the the plasma 
but its effects are macroscopic leading to a toroidal flow throughout the plasma column. 
Furthermore, according to Coppi's model for the tokamak case and our model for the mag- 
netotail, in both cases the cause is the presence of fluctuations due to pressure-gradient 
driven modes. 

However, we should warn the reader of some key differences that need to be considered in 
pursuing the analogy further. In geospace there is a toroidal magnetic field but it is small or 
comparable with the poloidal field. Furthermore, the plasma thickness, once rescaled with 
the ion inertial length, is smaller in the magnetotail (despite its physical size). The plasma 
beta is also different, being largely above unity in a large portion of the magnetotail. For 
these reasons the nature of the microinstabilities might be different and the analogy might 
not be complete. 

III. MODELING APPROACH 

To simulate the effect of small scale drift instabilities on the macroscopic momentum 
profile, we consider an idealization of the conditions typical of the current sheet in the Earth 
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magnetotail. The initial state is a Harris sheet pj3[, with magnetic field and plasma density 
given in GSM coordinates (x aligned with the main tail magnetic field in the Sun-Earth 
direction, y aligned with the current in the dawn-dusk direction, z aligned with the main 
gradients in the north-south direction) by B x = B t&iah(z/ L), n = n sech 2 (z / L) . 

The current sheet half-width is chosen as L = di/2, where di = c/u pi is the ion inertia 
length. The temperature is uniform and a temperature ratio Tj/T e = 5 and a mass ratio 
rrii/m e = 180 are assumed. To consider drift instabilities, we conduct the simulation in the 
2D (y, z) domain delimited by -6.4 < y/L < 6.4, < z/L < 25.6. 

To conduct the simulations, we rely on our long standing kinetic code CELESTE3D (see 
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Ref. |16j and references therein). The code is 3D but for the present study the x direction 
is not used. A grid with 128x256 cells is used with u P iAt = 0.1. In the y direction, periodic 
boundary conditions are used. In the z direction, the fields are subject to a Dirichlet bound- 
ary condition. Note that the latter condition allow the exchange of momentum between the 



plasma and the 
the solar wind) 



boundary (that summarizes the interaction with the Sun-Earth system via 



IV. SIMULATION RESULTS 

In the edge of the magnetotail current layer, the LHDI develops with properties described 



by previous theoretical studies [17|, |18|] . Its free energy is provided by the diamagnetic drift 
caused by the density gradient. The frequency of the LHDI is in the lower-hybrid range: 
Q ce > uj > Q C i. The LHDI is a micro instability with peak growth at scales of the order of 
the electron gyroradius: k±p e ~ 1 and propagating in the direction normal to the magnetic 
field: ku = 0. While the fastest modes are on the electron gyroradius scale and are mostly 

n □ 

electrostatic, longer wavelength (k± ~ y/p e Pi) electromagnetic modes exist [18|, |19j|. The 
LHDI remains mostly localized in the edge of the plasma because it is stabilized by the high 



beta in the centre of the layer 20] 



Figure [H shows the dawn-dusk (toroidal) current, electric field and velocity at time Vt c {t = 
14.7 when the LHDI is fully saturated. For each quantity, the left panel shows the 2D false 
color representation and the right panel reports the ?/-averaged profile (with the exception 
of the electric field fluctuation where instead of the mean, the RMS is shown). The presence 
of the LHDI is highlighted by the recurring stripes of red and blue in the toroidal electric 
field. The fluctuations peak at around one ion inertial length (corresponding to z = ±2L) 
off the center and remain localized in the edges, away from the center of the current sheet. 

The non-linear effect of the LHDI is to cause the spontaneous formation of a velocity 
profile in the dawn-dusk direction. The bottom panel of Fig. [1] shows the ion velocity at 
time Q C it = 14.7: the velocity presents a small scale noise due to the presence of the LHDI 
and even of Langmuir waves but on the whole it has developed a large scale structure absent 
in the original plasma. The shear is evident in the ^/-averaged profile, where the initial flat 
ion velocity profile is shown for reference. The saturation of the LHDI has caused a large 
exchange of momentum that has altered the profile producing a large shear. While the 
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FIG. 1: Simulation of the evolution of the LHDI, time shown: uj c it = 14.7. From top to bottom: 
toroidal current, electric field and ion flow speed. On the right, the quantities averaged in y are 
shown (solid blue line) compared with the initial profile (dashed red line). For the fluctuation 
toroidal electric field, the root mean square is shown. 



instability is limited to the edge, the momentum exchange does not remain localized in the 
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e layer, including its center. 
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221 ]: it is intensified and 



region where its cause is present but extends over the w 

The dawn-dusk current is also affected by the LHDI 
its initially straight surface is kinked. The kinking of the current layer has two separate 
scales with different causes. The smaller scale is still longer wavelength (k± ~ y/Pefk) than 
the bulk of the LHDI modes and it is due to the electromagnetic branch of the LHDI. 
Instead, the longer wavelength kinking is due to a secondary Kelvin-Helmholtz instability 
created by the formation of the velocity shear caused by the LHDI 111 ]. 

Evidence can be found in existing satellite databases to validate our results summa- 
rized above. As an example, Fig. [2] reports previously published results from the Geotail 
mission [231 ] . The time interval corresponds to one instance of the satellite entering the 
magnetotail current layer where the magnetic field reverses (at about time 11:07UT), visible 
in Fig. [2] when the magnetic field amplitude drops to near zero. As the satellite enters the 
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FIG. 2: GEOTAIL magnetic field data and plasma data provided by T. Nagai, H. Hayakawa and 
T. Mukai through DARTS at the Institute of Space and Astronautical Science (ISAS) in Japan. 
Data for March 30, 1995 when Geotail was at a distance of 15.2Re along the Earth magnetotail 



23]. 



layer, the plasma density increases, reaching its peak in the center and the temperature 
remains essentially constant within the layer. Therefore, the Harris equilibrium used in the 
simulations is a reasonable approximation of the real layer. The bottom panel shows the 
plasma speed, revealing the presence of a velocity shear with the peak-speed reached in the 
center and decreasing towards the edge. This finding is in remarkable agreement with the 
simulation results shown in the bottom panel of Fig. [TJ Additional Geotail data for the same 



23|, 



crossing shows a strong presence of electrostatic fluctuations in the edge of the plasma 
further supporting the scenario described above. 

A complete theory of the precise mechanism allowing the small scale LHDI developing in 
the edge to create a large velocity shear will require future investigations. However, one key 
aspect of it has been revealed in a recent simulation investigation 24J ]. highlighting the role 
of electrostatic fields. The formation of a velocity profile is accompanied by the formation 
of a large scale electrostatic field. Figure [3] shows for the same case in Fig. [1] the vertical 
component of the electric field. The saturation of the LHDI has created a large scale electric 
field that is primarily electrostatic in nature: the electric field shown is essentially curl-free 
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FIG. 3: Vertical electric field for the same run as in fig. [TJ time shown: iv c it = 14.7. On the left, 
the field is shown averaged in y. 



and it is caused by charge separation 24]. An obvious interpretation is that the small scale 
convective cells caused by the LHDI self-organize to form a large scale shear profile in a 



manner typical of many flow processes 



25). 



V. CONCLUSIONS 

We have reconsidered the recent theory by Coppi [9] relative to the analogy between the 
exchange of angular momentum in accretion disks and in tokamaks where a spontaneous 
toroidal rotation has been observed (e.g. Alcator C-Mod, DIII-D, TCV,Tore-Supra, JET, 
TFTR, COMPASS) . We have proposed above that a similar analogy can be 

extended to a new process recently discovered with a simulation investigation 10|, QjJ, [12] : 
the spontaneous formation of a velocity shear in the dawn-dusk direction in the Earth 
magnetotail by the action of drift waves (LHDI). We have outlined the possible analogy 
and we have presented a detailed simulation that investigates the conditions prevalent in 
the Earth environment. We have further showed satellite data that provide experimental 
evidence supporting the existence of the processes predicted by our simulations. 
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